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Abstract 

Two independent calculations of proton-deuteron elastic scattering observables including 
Coulomb repulsion between the two protons are compared in the proton lab energy region between 
3 MeV and 65 MeV. The hadron dynamics is based on the purely nucleonic charge-dependent AV18 
potential. Calculations are done both in coordinate space and momentum space. The coordinate- 
space calculations are based on a variational solution of the three-body Schrodinger equation using 
a correlated hyper spherical expansion for the wave function. The momentum-space calculations 
proceed via the solution of the Alt-Grassberger-Sandhas equation using the screened Coulomb 
potential and the renormalization approach. Both methods agree within 1% on all observables, 
showing the reliability of both numerical techniques in that energy domain. At energies below 
three-body breakup threshold the coordinate-space method remains favored whereas at energies 
higher than 65 MeV the momentum-space approach seems to be more efficient. 
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I. INTRODUCTION 


Although there is a long history of theoretical work on the solution of the Coulomb 


problem in three-particle scattering [l|, 


, the work of Refs 


• 3, E I pi 


pioneered the 


effort on fully converged numerical calculations for proton-deuteron (pd) elastic scattering 
including the Coulomb repulsion between protons together with realistic nuclear interactions. 
In their work the authors use the charge-dependent AV18 potential together with the Urbana 
IX three-nucleon force and proceed to solve the three-particle Schrodinger equation using the 
Kohn variational principle (KVP); the wave function satisfies appropriate Coulomb distorted 
asymptotic boundary conditions and is expanded at short distances in a pair correlated 
hyperspherical harmonics basis set. The results presented were fully converged vis-a-vis 
the size of the basis set and the angular momentum states included in the calculation. In 


parallel a benchmark calculation was performed j8[ where results obtained variationally were 
compared with those obtained from the solution of coordinate-space Faddeev equations for 
the AV14 potential at energies below three-body breakup threshold. 

In a recent publication p the momentum-space solution of the Alt-Grassberger-Sandhas 
(AGS) equation [l^ for two protons and a neutron was successfully applied, not only to 
pd elastic scattering but also to radiative pd capture and two-body electromagnetic disinte¬ 
gration of ^He. The treatment of the Coulomb interaction is based on the ideas proposed 


by Taylor 


Ilf for two charged particle scattering and extended in Ref. Q] for three-particle 


scattering with two charged particles alone. The Coulomb potential is screened, standard 


scattering theory for short-range potentials is used, and the o 
for the unscreened limit using the renormalization prescription 


in Ref. are converged vis-a-vis the screening radius R and the immber of included two 


otained results are corrected 
l[ [ll| . The results presented 


body and three-body angular momentum states. Although in Ref. the hadron dynamics 
is based on the purely nucleonic charge-dependent (CD) Bonn potential and its realistic 
coupled-channel extension CD Bonn -|- A, allowing for single virtual A-isobar excitation, 
other realistic potential models may be used easily as well. 

Motivated W recent experimental efforts in the measurements of pd elastic observ¬ 


ables 


nivatea ^ i 

HHQ, 


in the present paper we present benchmark results for a number of pd 


elastic scattering observables, both below and above three-body breakup threshold, using 


the charge-dependent AV18 m two-nucleon potential and no three-nucleon force. In Sec. E] 
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we make a short description of the methods we use, in Sec. IIIII we present the results, and 
in Sec. Cvl the conclusions. 


II. THE METHODS 


In this section we briefly introduce both methods and provide the basic framework for 
a general understanding of the technical procedures; further details may be found in the 
appropriate references. We choose to describe the method based on KVP using its traditional 
notation, which we attempt to carry over to the discussion of the integral equation approach 
in Sec. mu and uni Therefore the presentation of the integral equation approach will not 


be in the notation used in Ref. 




A. The Kohn variational principle 


The KVP can be used to describe nucleon-deuteron (Nd) elastic scattering. Below the 
three-body breakup threshold the collision matrix is unitary and the problem can be for¬ 
mulated in terms of the real reactance matrix (iP-matrix). Above the three-body breakup 
threshold the elastic part of the collision matrix is no longer unitary and the formula¬ 
tion in terms of the R-matrix, the complex form of the KVP, is convenient. Referring to 
Refs. Q,H,|i 7| for details, a brief description of the method is given below. The scattering 
wave function fw.f.') T is written as sum of two terms 


T = T 


c 




( 1 ) 


which carry the appropriate asymptotic boundary conditions. The hrst term, Tc, describes 
the system when the three-nucleons are close to each other. For large interparticle sepa¬ 
rations and energies below the three-body breakup threshold it goes to zero, whereas for 
higher energies it must reproduce a three outgoing particle state. It is written as a sum of 
three Faddeev-like amplitudes corresponding to the three cyclic permutations of the particle 
indices. Each amplitude yj), where Xj,yj are the Jacobi coordinates corresponding 

to the Tth permutation, has total angular momentum JJ^ and total isospin TT^ and is 
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decomposed into channels using LS coupling, namely 


Nc 

'hc(xi,yi) = '^(j)aixi,yi)yaijk,i) 

a=l 


ya{jk,i) 



( 2 ) 

(3) 


where Xi, yi are the moduli of the Jacobi coordinates and 3^^ is the angular-spin-isospin func¬ 
tion for each channel. The maximum number of channels considered in the expansion is Nc- 
The two-dimensional amplitude cfy^ is expanded in terms of the pair correlated hyperspherical 


harmonic basis 


la 


19| 


(l)a{Xi,yi) = p 




(4) 


. K 

where the hyperspherical variables, the hyperradius p and the hyperangle (/>*, are dehned 
by the relations Xi = p cos (pi and yt = p sin (pi. The factor ^‘^^P^^{(p) is a hyperspherical 
polynomial and fa{xi) is a pair correlation function introduced to accelerate the convergence 
of the expansion. For small values of the interparticle distance fa{xi) is regulated by the 
NN interaction whereas for large separations fa{xi) 1. 

The second term, in the variational wave function of Eq. o describes the asymptotic 
motion of a deuteron relative to the third nucleon. It can also be written as a sum of three 
amplitudes with the generic one having the form 


n 


LSj(Xi,y*)= Yl 


TT, 


(5) 


la=0,2 


where wi^{xi) is the deuteron w.f. radial component in the state la = 0,2. In addition, 
sP = 1, fP = 0 and L is the relative nucleon-deuteron angular momentum. The superscript 
A indicates the regular (A = R) or the irregular (A = I) solution. In the pd {nd) case, the 
functions TZ^ are related to the regular or irregular Coulomb (spherical Bessel) functions. 
The functions can be combined to form a general asymptotic state 


where 


^LSj(xi,yi) = f 2 °sj(xi,yi) PLU^l'S'A^i^yi) . 

US' 


( 6 ) 


^LSj(Xi,yi) = MoofJLSj(Xi,yi) +Moif^isj(Xi,yi) , 
^LSj(xi,yi)= Miof2f5j(xi,yi)+Miifi[gj(xi,yi) . 
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(7) 

( 8 ) 







The matrix elements Uij can be selected according to the four different choices of the matrix 
C = /^-matrix, i^“^-matrix, S'-matrix or T-matrix. A general three-nucleon scattering w.f. 
for an incident state with relative angular momentum L, spin S and total angular momentum 
J is 

'J'lsj = 5^ [^c(xi,yi)+ , (9) 

i=l,3 

and its complex conjugate is '^2sj- ^ variational estimate of the trial parameters in the 
w.f. can be obtained by requiring, in accordance with the generalized KVP, that the 

functional 

Yclf,] = ^Clf, - , (10) 

be stationary. Below the three-body breakup threshold, due to the unitarity of the S'- 
matrix, the four forms for the £-matrix are equivalent. However, it was shown that when 
the complex form of the principle is used, there is a considerable reduction of numerical 
instabilities j^. Above the three-body breakup threshold it is convenient to formulate the 
variational principle in terms of the S'-matrix. Accordingly, we get the following functional: 


= 'Sit' + - E\^ts'j) 


( 11 ) 


The variation of the functional with respect to the hyperradial functions u'^{p) leads to 
the following set of coupled equations: 


a' ,k' 


[^kk' (p)^ + ^kk' (p)^ + Cfcfc" (p) + ^2 


M, 


N 


E KS'ip) ui(fi) = 


( 12 ) 


For each asymptotic state two different inhomogeneous terms are constructed corre¬ 

sponding to the asymptotic ^Isj functions with A = 0,1. Accordingly, two sets of solutions 
are obtained and combined to minimize the functional m with respect to the S-matrix ele¬ 
ments. This is the hrst order solution, the second order estimate of the S'-matrix is obtained 
after replacing the hrst order solution in Eq. (HH). 

In order to solve the above system of equations appropriate boundary conditions must be 
specihed for the hyperradial functions. For energies below the three-body breakup threshold 
they go to zero when p —>■ cx), whereas for higher energy they asymptotically describe the 
breakup conhguration. The boundary conditions to be applied in this case have been dis¬ 
cussed in Refs. 0,ll3,E| and are briehy illustrated below. To simplify the notation let us 
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label the basis elements with the index n = [a,K], and introduce the completely antisym¬ 
metric correlated spin-isospin-hyperspherical basis element Q^{p,Q) as linear combinations 
of the products 

3 

, (13) 

i=l 

which depend on p through the correlation factor. In terms of the Q^{p-, fl) the internal part 
is written as 

Nm 

•tc = , (M) 

fl=l 

with Nm the total number of basis functions considered. The hyperradial functions Uf^{p) 
and ujf^{p) are related by an unitary transformation imposing that the “uncorrelated” ba¬ 
sis elements obtained by setting all the correlation functions /a(xj) = 1, form an 

orthogonal basis. Explicitly, the matrix elements of the norm N behave as 
r 

N„.{p)= da a^{p,n)>Q,.(p,n)^ n';^, + ^ + o(i/p’^) , torp^oo, (is) 

J P 

where, in particular, 

A'™' = j da s;(n)'s“.(n). (le) 

is diagonal with diagonal elements either 1 or 0. Therefore, some correlated elements 
have the property: Q^(p, fl) —> 0 as p —>■ oo. In the following we arrange the new basis in 
such a way that for values of the index p < Nm the eigenvalues of the norm are Affj, = 1 and 
for Nm + 1 < p < Nm they are A/)i = 0. 

For p ^ oo, neglecting terms going to zero faster than the asymptotic expression of 
the set of Eqs. m rotated using the unitary transformation dehned above, reduces to the 
form 




ft" / r- + 




+Oi^)\u;^ip) = 0 , (17) 


Mat ydp"^ p^ j , . ^ 

where E = and JC^ = G^j, + 3/2. Here is the grand-angular quantum number 

dehned as Gfj, = la + La + 2K and the matrix y is dehned as 

w ^ f da s“(sj)t X Qpa). (18) 


The dimensionless operator y originates from the Coulomb interaction as 

.. _ Mn 1 -|- Tj^z 1 + 

^ 2h?Q ^ cos0j 2 2 

i=l 


(19) 
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It should be noticed that = 0 if /i, /x' > Nm- 

In practice, the functions uj^{p) are chosen to be regular at the origin, i.e. a;^(0) = 0 and, 
in accordance with the equations to be satished for p —»• cx), to have the following behavior 
(h < Nm) _ 

Nm 

^ - E V e'*" . (20) 

fl’ = l 

where 6 ^/ are unknown coefficients. This form corresponds to the asymptotic behavior of 
three outgoing particles interacting through the Coulomb potential [2^. In the case of nd 
scattering (x = 0) the outgoing solutions evolve as outgoing Hankel functions H^^\Qp) 

K(p) ^ 

For values of the index p > Nm the eigenvalues of the norm are = 0 and the leading 
terms in Eq. (ED vanish. So, the asymptotic behavior of these functions is governed by 
the next order terms. However, for p > Nm, it is verified that ^ 0 as p —>• cx). 

In order to solve the system of Eqs. oa the hyperradial functions are expanded in terms 
of Laguerre polynomials plus an auxiliary function 

M 

E +‘IP„(P) , ( 21 ) 

m=0 

where z = jp and 7 is a nonlinear parameter. The linear parameters (m = 0,...., M + 1) 
are determined by the variational procedure. 

The inclusion of the auxiliary functions uJ^{p) dehned in Eq.(j2ID is useful for reproducing 
the oscillatory behavior shown by the hyperradial functions for p > 30 fm. Otherwise a 
rather large number M of polynomials should be included in the expansion. A convenient 
choice is to take them as the regular solutions of a one dimensional differential equation 
corresponding to the p-th equation of the system whose asymptotic behavior is the one of 
Eq-(ED- El the cases considered here the solutions obtained for the F-matrix stabilize for 
values of the matching radius po > 100 fm. 


B. The integral equation approach 


Q 


The integral equation to be solved is the AGS equation m for three-particle scattering 
where each pair of nucleons interacts through the strong potential v and the Coulomb 


potential w acts only between charged nucleons. The work in Ref. [9|] follows the seminal 
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work of Refs. [1|, lll| in the sense that the treatment of the Coulomb interaction is based 
on screening, followed by the use of standard scattering theory for short-range potentials 
and renormalization of the obtained results in order to correct for the unscreened limit. 
Nevertheless there are important differences relative to Ref. that are paramount to the 
fast convergence of the calculation in terms of screening radius R and the effective use of 
realistic interactions: 

a) We work with a screened Coulomb potential 


WR{r) = w{r) e hZ-R)'' 


( 22 ) 


where w{r) = j is the true Coulomb potential, a being the hne structure constant and 
n a power controlling the smoothness of the screeniM. We prefer to work with a sharper 
screening than the Yukawa screening (n = 1) of Ref. p because we want to ensure that the 
screened Coulomb potential wr approximates well the true Coulomb one w for distances 
r < R and simultaneously vanishes rapidly for r > R, providing a comparatively rapid 
convergence of the partial wave expansion. In contrast, the sharp cutoff {n —> oo) yields 
unpleasant oscillatory behavior in momentum space representation, leading to convergence 
problems. We hnd values 3 < n < 6 to provide a sufficient smoothness and fast convergence; 
n = 4 is used for the calculations of this paper. 

b) Although the choice of the screened potential improves the partial wave convergence, 
the practical implementation of the solution of AGS equation still places a technical difficulty, 
i.e., the calculation of the AGS operators for nuclear plus screened Goulomb potentials 
requires two-nucleon partial waves with pair orbital angular momentum considerably higher 
than required for the hadronic potential alone. In this context the perturbation theory 
for higher two-nucleon partial waves developed in Ref. Q is a very efficient and reliable 
technical tool for treating the screened Goulomb interaction in high partial waves. 


As a result of these two technical implementations, the method 


24l | that was developed 


before for solving three-particle AGS equations without Goulomb could be successfully used 
in the presence of screened Goulomb. Using the usual three-body notation, the full multi¬ 
channel transition matrix reads 

<’(Z) = VTy>(Z)G„(Z)C/™(Z), (23a) 

cr 

where the superscript {R) denotes the dependence on the screening radius R of the Goulomb 
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potential, Go(Z) = (Z — Hq) ^ the free resolvent, = 1 — ^/ 3 a; and 

Ti^\Z) = (n„ + w^r) + (n„ + Wo,R)G,{Z)Tj^\Z). (23b) 

The two-particle transition matrix T^\z) results from the nuclear interaction Va between 
hadrons plus the screened Coulomb WaR between charged nucleons {waR = 0 otherwise). As 
expected the full multichannel transition matrix (Z) must contain the pure Coulomb 
transition matrix T^^{Z) derived from the screened Coulomb between the spectator 

proton and the center of mass (c.m.) of the remaining neutron-proton {np) pair in channel 
a 


RS(Z) = + »UTG™(Z)T„X(Z). 


(24) 


where — 0 for n{pp) a channels and the channel resolvent 

G™(Z) = (Z-Ha- - w^n)-K 


(26) 


In a system of two charged particles and a neutral one, when WaR = 0, W^r ^ 0 and vice 
versa. 


As demonstrated in Refs. |1|, |2| the split of the multichannel transition matrix 


<>(Z) = aa„T‘S‘(Z) + V'SPZ) - '5a=T„T(Z)| 


(R)/ 


(26) 


into a long-range part 6 j 3 aT^R{Z) and a Coulomb distorted short-range part \U^^{Z) 
^/ 3 a^ai?^(^l_is extremely convenient to recover the unscreened Coulomb limit. According 


m 


to Refs, [ll, |2| the full pd transition amplitude {(j)p{(lf)i^Pf\Ui 3 a\(j>ai^i)i^ai) for the initial and 
hnal channel states with relative pd momentum q* and q/, qj = qi, energy Ea{qi), and 
discrete quantum numbers iZq,. and z/^^, is obtained via the renormalization of the on-shell 
U^^{Z) with Z = Ea{qi) -l- iO in the inhnite R limit. For the screened Coulomb transition 
matrix T^^{Z), contained in UpJ{Z), that limit can be carried out analytically, yielding the 
proper Coulomb transition amplitude (0^(q/)z/;3j.|T^^|0Q,(qj)z/Q,J [ll,llli, while the Coulomb 
distorted short-range part requires the explicit use of a renormalization factor, 

( 0/3 (q/)I 10 a ) 

= <5/3a(0/3(q/)z^/3^|T^C I0a(q*)i^a0 

+ Jta/^{ZR“(9/)('3u(q/)>'/J,|[C/9f’(£^o(9i) + «) 

- + *0)l|,/.„(q.)R„,)Z«b9.)}. 


( 27 ) 
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The renormalization factor 


ZR{q) = e 




(28a) 


contains a phase (pniq) which, though independent of the pd relative orbital momentum L 


in the inhnite R limit, is given by 
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= (^dq) - VLR{q), (28b) 

where rjLRiq) is the diverging screened Coulomb phase shift corresponding to standard 
boundary conditions, and o'L(g) the proper Coulomb phase referring to logarithmically dis¬ 
torted Coulomb boundary conditions. The limit of the Coulomb distorted short-range part 
of the multichannel transition matrix — 6 fSaTa^{Z)] has to be performed numeri¬ 

cally but, due to its short-range nature, the limit is reached with sufficient accuracy at hnite 
screening radii R. Furthermore, due to the choice of screening and perturbation technique 
to deal with high angular momentum states, (Z) — 6 j 3 aT^R{Z)] is calculated through 
the numerical solution of Eqs. (121 and ( 121 , using partial-wave expansion. 

In actual calculations we use the isospin formulation and, therefore, the nucleons are 
considered identical. Instead of Eq. (1^ we solve a symmetrized AGS equation 

f/(-R)(Z) = PGq\Z) + PT^^\Z)Go{Z)U^^\Z), (29) 


P being the sum of the two cyclic three-particle permutation operators, and use a properly 
symmetrized pd transition amplitude 


(0a(q/)I f/|0a(q*)J 

= (0a(q/)^2a^|(Cc I0a(q0i^ai) 

+ Jim {^R"(g/)(0a(q/)z2a^|[f/^-^H^a(?0 +^0) 

R—>oo 

for the calculation of observables. For further technical details we refer to Ref. 


(30) 


III. RESULTS 


In this section we compare numerical calculations for a number of elastic observables 
performed using the KVP and the integral equation approach. Three different lab energies 
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have been considered: 3, 10, and 65 MeV. The Coulomb effects are expected to be sizable 
in most of the observables at the first two energies. The two methods use a different scheme 
to construct the scattering states with total angular momentum and parity In the 
KVP the LS coupling is used and channels are ordered by increasing values of £„ + Lq,. 
The expansion of the scattering state is truncated at values ia + La = L^ax + 2, where 
Lmax is the maximum value of L corresponding to the asymptotic states In the 

integral equation approach the jj scheme has been used. The channels have been ordered 
for increasing values of the two-body angular momentum j and for the strong interaction the 
maximum value jmax = 5 has been considered for the first two energies whereas at 65 MeV 
the value jmax = 6 has been used; the screened Coulomb interaction is taken into account 
up to jmax = 25 as described in Ref. j^. Both numerical calculations presented here are 
converged relative to the number of three-body partial waves. In addition, the variational 
calculations are converged relative to the size of the hyperspherical basis set and, in the 
integral equation approach, the results are converged with respect to the screening radius 
R. 


In Figs. Hand El we compare the differential cross section and vector and tensor analyzing 
powers for pd elastic scattering at the three selected energies, 3, 10 and 65 MeV proton 
lab energies. In Fig. El a selection of spin transfer coefficients at 65 MeV is shown. In the 
figures, two different curves are shown corresponding to calculations using the KVP (thin 
solid line) and integral equation approach (dotted line). By inspection of the figures one may 
conclude that the agreement is excellent because the numerical calculations agree to better 
than 1%. In fact the curves are practically one on top of the other, the exceptions being the 
maximum of T 21 and some spin transfer coefficients at 65 MeV in which a small disagreement 
is observed. Nevertheless, it is important to mention that in all cases the difference between 
the two curves is smaller than the experimental accuracy for the corresponding data sets. 
Likewise the agreement between the two calculations largely exceeds the agreement of any 
of them with the data as shown in Refs. jJ, 

The present results can be used to study Coulomb effects by comparing nd to pd calcu¬ 
lations. In Fig. |3]we analyze the evolution of the Coulomb effects for the differential cross 
section, the nucleon analyzing power Ay and two tensor analyzing powers, T 20 and T 21 at 3, 
10 and 65 MeV proton lab energies. In order to reduce the number of curves in the figure 
for the sake of clarity we present results obtained using the integral equation approach. The 
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results obtained using the KVP for nd scattering agree at the same level already shown for 
the pd case in the previous hgures. In Fig. |3]the thin solid line denotes the pd calculation 
whereas the dotted line denotes the corresponding nd calculation. The latter agrees well 
with the results of other existing nd calculations j^. From the hgure we observe that 
Coulomb effects are appreciable at 3 and 10 MeV but are considerably reduced at 65 MeV. 


A more exhaustive analysis on Coulomb effects can be found in Refs, |9|, |2( 

In addition to the benchmark comparison using AV18 potential we also give one result for 
the Malfliet-Tjon (MT) 1-111 potential, in order to resolve an existing problem. Reference 


Q 


reports a disagreement between pd phase shifts results for MT 1-111 potential calculated 


using the hrst technique of this paper, the KVP |^, and the conhguration-space Faddeev 


equations [ 3 ]. The calcuMion based on the second technique of tMs paper, the momentum- 


Q. 


space integral equations j^, clearly confirms the results of Ref. A detailed comparison 
of pd and nd phase shift results for MT 1-111 potential is given in Table U 

In the following we discuss some of the limitations inherent to the two methods used to 
describe pd elastic scattering. The KVP, as presented here, reduces the scattering problem to 
the solution of a linear set of equations in which the matrix elements of the Hamiltonian have 
to be computed between basis states; increasing the energy, appreciable contributions from 
states with high values of + La appear. In order to take into account these contributions, 
a very large basis has to be used with the consequence that numerical instabilities start 
to appear. In the integral equation approach at very low energies convergence in terms of 
screening radius requires R > 30 fm, which in turn increases the number of two-body partial 
waves that are needed for convergence. The interplay of these two requirements makes 
the integral equation solution unstable at those very low energies. An interesting heuristic 
argument to understand the size of the screening radius needed for convergence is the wave 
length A corresponding to the on-shell momentum. At 3 MeV, 10 MeV, and 65 MeV proton 
lab energy, for which a screening radius of 20 fm, 10 fm, and 7 fm is needed for convergence, 
A is 24.8 fm, 13.6 fm, and 5.3 fm, respectively. It appears that for the calculation of pd 
elastic scattering observables the screening has to be only so large that one wave length can 
be accommodated in the Coulomb tail outside the range of the hadronic interaction; seeing 
proper Coulomb over one wave length appears enough to provide, with the additional help 
of renormalization, the true Coulomb characteristics of scattering despite screening. 


Q. 
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IV. CONCLUSIONS 


In the present paper two methods devised to describe elastic pd scattering are compared 
for a wide range of energies. One of the methods, the KVP, was developed a few years 
ago and used to study how realistic potential models, including two-body and three-body 
forces, describe the elastic observables measured for that reaction. On the other hand, 
numerical accurate results have been recently obtained solving the AGS equation for pd 
scattering using a screened Coulomb potential corrected for the unscreened limit using a 
renormalization prescription. As has been briefly described in the present paper, both 
methods are substantially different. It is satisfactory to observe that both methods produce 
essentially the same results for a large variety of elastic observables using a realistic two- 
nucleon potential. We stress the fact that the selection of observables here presented is only 
part of the observables compared. In all cases, similar patterns have been obtained. In 
addition, by comparing the pd calculations to the corresponding nd calculations. Coulomb 
effects have been estimated. As expected these effects are sizable at low energies but at the 
highest analyzed energy, 65 MeV, they are small, except at forward scattering angles. From 
these considerations it is possible to identify on a hrm basis which pd observables may or 
may not be analyzed by calculations in which the Coulomb interaction has been neglected. 

We can conclude that at present it is possible to describe pd elastic scattering, including 
the Coulomb repulsion, using standard techniques as the Faddeev equations in conhguration 
and momentum space or variational principles. Moreover, in Ref. j27| the treatment of other 
terms of the NN electromagnetic potential as the magnetic moment interaction has been 
discussed. 
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to 

o 

^Vo 


4 

^0 

nd at 14.1 MeV 

105.48 

0.4649 

68.95 

0.9782 


105.49 

0.4649 

68.95 

0.9782 

nd at 42.0 MeV 

41.34 

0.5022 

37.72 

0.9033 


41.36 

0.5022 

37.71 

0.9033 

pd at 14.1 MeV 

108.44 

0.4984 

72.60 

0.9795 


108.39 

0.4983 

72.62 

0.9796 

pd at 42.0 MeV 

43.67 

0.5056 

39.95 

0.9046 


43.70 

0.5056 

39.97 

0.9046 


TABLE I: nd and pd phase-shift and inelasticity parameters calculated with MT I-III potential. 
For each reaction, the KVP and integral equation approach results are given in the first and second 
row, respectively. 
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FIG. 1; Differentia[ cross section, and the proton and deuteron anatyzing powers Ay and iTn are 
shown for pd elastic scattering at 3 MeV (first row), 10 MeV (second row) and 65 MeV (third row) 
lab energies as functions of the c.m. scattering angle. Results obtained using KVP (thin solid line) 
and integral equation approach (dotted line) are compared. 
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FIG. 2; Tensor analyzing powers T 20 , T 21 , and T 22 are shown for pd elastic scattering at 3 MeV 
(first row), 10 MeV (second row) and 65 MeV (third row) lab energies as functions of the c.m. 
scattering angle. Results obtained using KVP (thin solid line) and integral equation approach 
(dotted line) are compared. 
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FIG. 3: A selection of proton to proton and proton to deuteron spin transfer coefficients for pd 
elastic scattering at 65 MeV lab energy as function of the c.m. scattering angle. Results obtained 
using KVP (thin solid line) and integral equation approach (dotted line) are compared. 
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FIG. 4; Calculations for pd (solid line) and nd (dotted line) scattering are compared for differential 
cross section, Ay, T20 and T21 at 3 MeV (first column), 10 MeV (middle column) and 65 MeV (last 
column) nucleon lab energies. 






